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Abstract. The Casimir forces between two plates moving parallel to each other at arbitrary 
constant speed are found by calculating the vacuum electromagnetic stress tensor. The 
perpendicular force between the plates is modified by the motion but there is no lateral force 
on the plates. Electromagnetic vacuum fluctuations do not therefore give rise to "quantum 
friction" in this case, contrary to previous assertions. The result shows that the Casimir-Polder 
force on a particle moving at constant speed parallel to a plate also has no lateral component. 

PACS numbers: 12.20.Ds, 42.50.Lc, 46.55.+d 

1. Introduction 

Over half a century ago, the vacuum force between two parallel dielectric plates was 
calculated by Lifshitz |fTl!2l. Casimir's simple formula [3| for the vacuum force between 
two perfect mirrors then emerges in the limit of infinite permittivity. A generalization of 
the problem to where one of the plates moves at a constant speed parallel to the other 
might be expected to be straightforward, but in fact authors have obtained conflicting 
answers for the vacuum forces in this case [HI [51 |6l d [H HI \W[ [HI [121. Interest in this 
problem has attached particularly to the lateral component of the Casimir force, which 
would give information on the extent of a quantum vacuum contribution to friction between 
bodies [^|5l|6l[7l[8l[9l[IOl[IIl[l2l[l3l[l4l[Il. In this paper we use Lifshitz theory [|Il[l6l[l to 
find the zero-temperature Casimir forces between the plates for arbitrary velocity. Despite a 
general acceptance in the literature on this problem [4l[5l|6l[7l|8][9l[T0l[ITl[T2|| that it provides 
an example of quantum friction, the exact solution, presented here for the first time, shows 
that the motion does not in fact induce a lateral Casimir force. 

The geometry and notation of the problem are set out in Figure[T]and caption. The plates 
have arbitrary permittivities and permeabilities and are taken to have a separation a in the x- 
direction. Plate 2 moves in the positive y-direction and we use /? to denote its speed in units of 
the speed of light c, i.e. the speed is /3c. It is clear that the arrangement in Figure [T] represents 
one of the most obvious and basic problems in the theory of Casimir forces and the relatively 
limited attention it has received is somewhat surprising (a dozen papers or so, over the past 
thirty years). The practical significance of this problem has also increased with the emerging 
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Figure 1. Two plates with general electric permitivities and magnetic permeabilities lie in the 
2/z-plane with constant separation a. The second plate moves in the postive y-direction with 
speed /3c, where c is the speed of light. We compute the Casimir forces on Plate 1 . 

importance of Casimir forces in micro- and nano-engineering ifTTlfTSll . Indeed, the problem 
in question is the simplest possible model of a nanomachine. 

Before entering into technical details it is worth considering what one might reasonably 
expect for the Casimir forces in this case. Casimir forces are caused by the vacuum zero-point 
modes of the electromagnetic field [19]. If the plates in Figure [T] are perfect mirrors then 
the zero-point modes do not penetrate the materials and the motion of Plate 2 can therefore 
have no effect on the vacuum forces: there is an attractive force between the mirrors given by 
Casimir's formula [|3l with no lateral component. For realistic materials, however, the zero- 
point modes will penetrate the moving plate and the motion should therefore affect the Casimir 
force; the perpendicular force is indeed altered by the motion, but the intriguing question is 
whether the vacuum modes are influenced in a way that produces a lateral force on the plates. 
At first glance a lateral force may seem plausible: the quantum vacuum can sometimes be 
usefully thought of as a kind of fluid [15J, and one could imagine this fluid exerting a drag 
on the moving plate, with the reciprocal drag of the vacuum exerting a lateral force on the 
non-moving plate. The quantum vacuum is a very strange kind of fluid, however: uniform 
motion relative to the quantum vacuum has no meaning because of the Lorentz invariance 
of quantum field theory, so an isolated moving body only experiences a quantum-vacuum 
"viscosity" when it accelerates, not when it moves at a constant velocity [|20l . Whether the 
presence of the non-moving plate in Figure 1 is enough to make the quantum vacuum become 
viscous to the moving plate (and to the non-moving one) is far from clear. The following 
consideration suggests the absence of a lateral force. A moving medium is equivalent to 
a particular non-moving bi-anisotropic medium ||2T1 1221 : for the present argument the bi- 
anisotropy can be very small, so it is hard to see any reason why such a material could not be 
constructed in principle. It would be very strange if a bi-anisotropic medium could be used 
to induce a unidirectional lateral Casimir force as this would seem to allow the extraction of 
unlimited energy from the quantum vacuum. If anisotropy is introduced into both plates then 
lateral Casimir forces are certainly possible, but these are rather torques that act to orient the 
plates to an equilibrium position Il23ll24ll25ll26l . 

Almost all of the treatments of this problem [il[51[6l|7l[8l|9l[l0l[IIl[l2l have made 
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use of various approximations, and there is disagreement on the magnitude of the Casimir 
forces. We make particular mention of the paper by Barton [|9l, although the concern in 
that analysis was almost exclusively with the case where the moving plate in Figure 1 is 
accelerating rather than moving at a constant speed. The approximations used in f9l included 
ignoring the dispersion and absorption of the plates and considering only the lowest-order 
effects of the motion. The lateral force on the plates was found to vanish in the case of constant 
velocity. Barton noted that this result contradicted previous treatments of the constant- 
velocity case in BH m m (which also contradict each other) but because of the different and 
severely limited methods used to model the materials he drew no conclusions as to the correct 
result for realistic materials. All other treatments of the problem with which we are concerned 
conclude that there is a lateral force on the plates. The previous claims 111211 (see also [l35l ) of 
an exact solution were based on Rytov's theory of electromagnetic fluctuations, which utilizes 
the classical Green tensor for the problem to calculate the forces. Since our results are also 
derived from the classical Green tensor it is here that our work makes direct contact with the 
existing literature. As we show in Section |4} the Green tensor was not calculated correctly 
in [lT2l . even approximately (essentially the same error was also made in llTOl ). In fact, the 
classical Green tensor represents the main computational challenge in this problem and it is 
presented here for the first time. 

In Section [2] we show how the vacuum electromagnetic stress tensor is related to the 
classical Green tensor; this is similar to the usual Lifshitz theory ^ [161 13 except that we 
treat only zero temperature. The exact Green tensor for the problem is derived in Section [3| 
this is achieved by a novel method which utilizes physical reasoning in place of considerable 
algebraic labour. Section |4] presents the Casimir force on the plates, and in our Conclusions 
we discuss the implications for the Casimir-Polder force between a particle and a plate. 



2. Lifshitz theory at zero temperature 

Casimir forces are given by the vacuum expectation value of the electromagnetic stress 
tensor 023: 

a- = £o(E ® E) + /io ^(B ® B) - ^l(£o(^') + /io '(^'))- (1) 

The methodology of Lifshitz theory is to compute the expectation values in the vacuum stress 
([T]) using the retarded Green tensor of the classical vector potential A in a gauge in which the 
scalar potential set to zero 0. In this gauge the electric and magnetic fields are given by 

E = -dtA, B = V X A, (2) 

and the vector-potential wave equation is 

Vx Vx+i^^^ A = /ioj, (3) 

where j is the current density of any sources . The retarded Green tensor for the vector potential 
satisfies 

(^V X V X +^5') G(r, t- r', t') = 16{r - r')6{t - t'), (4) 
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and relates the vector potential at r, t to the current density j at r', t' , where t' < t: 

A(r,t)=/io r dt' I dVG(r,t;r',t')-j(r,t')- (5) 

J — oo J 

The retarded boundary condition means that G(r,t; r',t') = for t' > t and in frequency 
space the Green tensor is defined by 

POD 

G(r,r',u;)=/ dt G(r, t, r', 0) e^^*, (6) 
Jo 

VxVx — - G(r,r',cj) = 15(r-r'). (7) 



c 



Since a real current will produce a real vector potential, ([5]) shows the Green tensor is real and 
in frequency space this gives the property 

G(r,r',-^*) = G*(r,r',^). (8) 

The similar property e{—lu*) = e*{uj) of the dielectric permittivity is derived in an identical 
manner [|27]| . 

In (|4]) we have written the equation for the Green tensor in the vacuum between the plates; 
the Green tensor inside the plates will not be required to impose the boundary conditions on 
(|7]) since we reduce the problem to a consideration of plane waves, the boundary conditions for 
which are imposed by their standard reflection and transmission coefficients fi27il . Equation 
([T]) is also valid only in the vacuum between the plates. The net force on the plates is a result 
of the electromagnetic stress tensor in the plates as well in the vacuum but it turns out that the 
relevant stresses in the plates are zero, as is discussed in detail in Section |4] and Appendix B 



At finite temperature it can be shown through use of the fluctuation-dissipation 
theorem and an assumption of time- symmetric boundary conditions, that the retarded Green 
tensor has a simple relation to a correlation function for the vector potential. Clearly this 
derivation does not apply to the problem considered here. For zero-temperature problems 
there is in fact no need to invoke the fluctuation-dissipation theorem and there is a much 
simpler proof [28 1 of the required relations. 

The vector potential operator A(r, t) is given by ||29ll 



A(r,t) = [Ak,.(r)e-^-'=*ak,. + A*,,(r)e 

k,o- 



(9) 



where Ak,cr(r), A^^(r) are a complete set of modes and a labels two linearly independent 
polarizations. We use the standard technique [29] of a finite quantization cavity, the volume of 
which can eventually be increased without hmit. The modes are solutions of the homogeneous 
equation 



2 " 

V X V X -^^ 1 Ak,.(r) = 



c2 



and have the normalization 



A^,,(r) ■ Ak',a'(r) d^a; = ^^l^^fc"^--'' / Ak,a(r) • Ak',a'(r) d'x = 0. (10) 
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Figure 2. Contours in the complex w-plane. The left picture shows the contour C. In the 
centre picture the integration below the negative real axis has been replaced by an additional 
integration above the positive real axis, as described in the text. In the right picture the contour 
above the poles on the positive real axis has been rotated to the imaginary axis. 



The Green tensor (|6])-(|7]) can be expanded in terms of the eigenfunctions (|9]) by a standard 
procedure (see [|27ll for example). A (positive frequency) Green tensor G{y,v' ,uj) with the 
same boundary conditions as the modes Ak,cr(i") can be expanded in terms of these modes. 
We can therefore write 

G(r,r',^) = ^ Ak,.(r) ® Ck,,(r'). (11) 

k,o- 

Acting on ( [TT| ) with the operator on the left-hand side of (|7]) we obtain 
5^ Ak,. ® Ck,a = 15(r - r'). 



A dot product with A^^(r') and use of (10) then determine Ck,cr and we find the standard 
expansion [i27ll 

G(r, r , ^) = ^ E ^^Ak,.(r) ® A^A^'). (12) 

Consider the equal-time vacuum correlation function for the electric field: (E(r, t) ® 
E(r', t)). Using ^ and the operator version of ^ we find 

(E(r, t) ® E(r', t)) =J2^l Ak,.(r) ® A*_,(r'). (13) 



The right-hand side of ( 13 1 is proportional to a contour integration of the Green tensor ( 12 1 
in the complex u plane. The denominator cu^ — cu^ = {cuk — to) (tOk + to) gives poles in the 
complex CO plane at all cok and —10^. Consider a contour C that lies just below the poles at —to^ 
on the negative real axis and just above the poles —cok on the positive real axis and which is 
then closed in the upper half complex plane (see Figure|2]). From ([6]), G(r, r', to) is analytic 
in the upper half plane and vanishes for Im(a;) 00. An integration of ( [T2] ) over the contour 
C therefore has contributions from the —cOk poles only. In this manner we find 

(E(r, t) ® E(r', t)) = / G(r, r', ^)a;2d^. (14) 

In order to compute the stress tensor ([T]) we require the correlation function ( 14) at r = r'. The 
electric field operators commute with each other at space-like separated events, but the limit 
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r ^ r' is singular and must be regularized (see below). To avoid ambiguities in the operator 
ordering we follow the usual prescription of quantum mechanics and define the correlation 
function at r = r' by 



(Ei(r,t)®4-(r,t)) 
We then obtain from ([14]) 

(E(r,t)®E(r,t)) = 



lim - 

r^r' 2 



h 



{Ei{r,t)®Ej{r',t)) + {E,{r' ,t) ® Ei{r,t)) 



2e^& ivr Jc 



G fr, r, lS)ijJ dcj, 



where G'^ is the symmetrized Green tensor, given by 



1 



r,r 



Gji(r',r,u;)] 



(15) 



(16) 



Since G(r, r',u;) vanishes for Im(cj) oo (see ([6])), the upper part of the semicircle 
in Figure [2] does not contribute to the integral (15). Moreover the contribution from large 
real u) does not contribute either due to the regularization procedure required for Casimir 
calculations in which the bare Green tensor in the absence of the plates is subtracted from 
G(r, r', a;) [l2l; this subtraction gives zero at large frequencies where the plates are transparent 
and macroscopic electromagnetism breaks down [,2J. Hence there is no contribution from the 
entire semicircle in Figure |2] 

Using ([8]), the integration in ( [T5| ) over the part of the contour C that lies along the negative 
real cj-axis can be replaced by an integration along the positive real cj-axis (see Figure [2]). 
Since the contour lies above the poles on the positive real axis, it can be pushed up to the 
imaginary axis (see Figure |2]). Eqn. ([8]) shows that the Green tensor is real for imaginary uj 
and we are left with the following integration over positive imaginary frequencies uj = \^: 



(E(r,t)®E(r,t)) 



TC 



G^(r,r,iOrde 



(17) 



The equal-time correlation function for the magnetic field operator is obtained in a similar 
manner using (|2]) and (|9]), with the result 



(B(r,t)®B(r,t)) 



TT 



lim 

r— >r' 



V X G*(r,r',iOx V d^, 



(18) 



where x V' denotes a curl on the second index of the Green tensor, so that for a vector V(r') 

we have Vx V'= V x V. The fact that the integrations in ( 17 ) and ( 18 ) are over imaginary 
frequncies (u; = iQ ensures that the integrals are well behaved LLi2iLl6J. The stress tensor 
([T]) follows immediately from ( 17 ) and ( 18 ). 



The expressions ( 17 ) and ( 18 ) are familiar from Lifshitz theory [|Tl [161 13, except that 
here the symmetrized Green tensor G"^ appears in these correlation functions rather than G. 
In the static case of non-moving plates there is in fact no difference since the Green tensor is 
diagonal when r = r' (in the Cartesian coordinates of Figure [T]). For the problem considered 
here, however, the emergence of the symmetrized Green tensor in ( [17] ) and ([TSj) ensures that 
the electromagnetic stress tensor is symmetric. We briefly return to this point in Section |4} 

Because of the interaction of the modes with the plates, which are in general dissipative. 



one may object that they cannot therefore be normalized as in ( 10). This technical issue of 
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quantum electrodynamics in absorptive media can be dealt with using Schwinger's source 



theory and its modern version OOl : here we used the modes as a tool to obtain (17) and 
( [T8| ), which remain valid in absorptive media ll30l . The issue of dissipation has also been 
much discussed in the Casimir literature (see and references therein) with the result that 
expressions for the Casimir force in terms of reflection and transmission coefficients (such 
as are derived here) are valid for dissipative media. This is intuitively reasonable since the 
quantum vacuum persists even in absorptive materials. 



3. The Green tensor 



With retarded boundary conditions, the solution of the monochromatic equation (|7]) has a 
simple physical meaning: an oscillating dipole at the point r' emits electromagnetic waves 
of frequency uj and G(r, r', uj, (3) is the resulting vector potential at the point r. The second 
index in Gij represents the orientation of the dipole at r', while the first index represents 
the components of the vector potential at r. Using this physical consideration it is clear from 
Fig.[T]that the solution will be a linear superposition of waves that have reflected off the plates, 
with the number of reflections ranging from zero to infinity. To write down the solution we 
exploit the homogeneity of the problem in the y- and z- directions to Fourier transform the 
Green tensor as follows: 



G(x, x', M, f , i^) 



(\y / d2;G(r,r',iOe 



-\u(y-y')-iv(z-z') 



(19) 



so that we decompose the waves emitted by the dipole into plane waves. In the absence of the 
plates the solution is the bare Green tensor ^ 



c 



w 
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^2 _|_ ^2 _|_ ^2 



(20) 



(21) 



(22) 



The two possibilities in (20) are plane waves propagating to the right (first line) or to the 
left (second line), with wave vectors k = {±.iw^u^v). The imaginary x-component of the 



wave vectors is a consequence of the imaginary frequency, and in (22) we simply have the 
relation u = ck. In physical terms the vacuum solution ([20]) is trivial: it is the only way 



the dipole can propagate plane waves from x' to x. In the presence of the plates both plane 



waves in (20) will reflect off the plates and reverse direction, so the left-moving plane wave 
can propagate from x' to x even if x > x', with similar considerations applying to the right- 
moving plane wave. Consequently, both the right- and left-moving waves will appear in the 
solution regardless of whether x is greater or less than x', in contrast to the vacuum solution 



(20). Let R2 be the reflection operator (matrix) that transforms a right- moving plane wave 
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at plate 2 into the resulting reflected left-moving plane wave, and let Ri be the reflection 
operator that transforms a left-moving plane wave at plate 1 into the reflected right-moving 
plane wave. We can now write down the solution (functional dependences are suppressed): 



G 



(1 
(1 



-2wa 



R1R2) 



(23) 



The first line in ( [23] ) subtracts the left- or right- moving plane wave in ( [20] ), depending 
on whether x is greater or less than x' . This subtraction is necessary because the direct 
propagation, without reflections, of both the right- and left-moving plane waves from x' to 



X is contained in the remaining terms in (23 1; but only one of these propagations is possible. 



depending on whether x is greater or less than x' , and the first line in ( [23] ) automatically 
subtracts the irrelevant one. The inverse matrices in ( [23] ) are geometric series representing 
every possible number of double reflections off both plates, the exponentials providing the 
propagation distance 2a for each double reflection. The initial right- and left-moving plane 
waves that leave x' reach x after both an even and odd number of reflections; this explains the 



terms multiplying the inverse matrices in (23). Each term in (23), after the series expansion 
of the inverse matrices, has an overall exponential factor that accounts for the propagation 
distance involved, with e^""*, s > 0, representing a propagation distance s to the right for the 
initial right-moving plane wave, but to the left for the initial left-moving plane wave. Note 



that we have not omitted any part of the derivation of ( [23] ); it was written down using the 
reasoning described above. 

It remains to deduce the reflection operators Ri and R2. Reflection of left-moving plane 
waves at plate 1 can be calculated by decomposing the plane wave into two components, one 
with polarization in the plane of incidence (TM) and the other with polarization perpendicular 
to this plane (TE); these components must then be multiplied by the well-known reflection 
coefficients for these polarizations ll27l . As well as the overall change in amplitude of the 
two polarization components, the x-component of A (proportional to E) changes sign on 
reflection, to remain perpendicular to k. The reflection operator Ri is therefore 



Ri 



/ -1 \ 

10, 

1 y 



v 



/il(lCKjW — Wi 



TBI 



ei{icK,jw — w\ 



Hi [ICK )W + Wi El {iCK)w + Wi 

Wi = a/m^ + t>2 + ei(ic/t)/ii(icK)fi;2. 



(24) 
(25) 

(26) 
(27) 



In (24), P^i (Pei) projects E to its component in (perpendicular to) the plane of incidence, 
while Hx flips the x-component of E (the i?-polarization has no x-component). Eqns. ( [26] ) 
are the reflection coefficients [|27ll . The (left- moving) wave vectors in the vacuum (k) and the 
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plate (ki) are 

k = (— iw, M, t>), ki = {—iwi,u,v). 
The projection operators are given by 

PEi = nEi^nEi, PBi = nBi®nBi 

n^i 



\ u J 



K\Jv? + 



( 



—uw 
—vw 



(28) 
(29) 
(30) 



where n^;! and n^i are unit vectors in the two polarization directions, for a left-moving plane 
wave. 

If /3 = 0, the reflection operator R2 for right-moving plane waves at plate 2 is similar to 



(p4))-([30|). As the waves are now right- moving, instead of ( [28] ) we have the wave vectors 

k = {iw, u, v), k2 = {1W2, u, v). 



We therefore simply change the material quantities in (24)-(30) to those of plate 2 and flip 



E' 



the sign of the w's (this affects the polarization direction pO|)). When /3 7^ this result for 
plate 2 holds in an inertial frame co-moving with the plate. To obtain R2 in the frame of 
Fig. [T] we must Lorentz boost in the negative y-direction. Denoting by primes any quantities 
in the co-moving frame that differ from their values in the frame of Fig. [TJ we have [|27ll 

-f{K + i(3u), u' = -f{u-i(3K), 7 = (l-/?2)-^^ (31) 

= ^{E^ + (3cB,), E', = ^{E,-(3cB,), (32) 

5; = 7(5, - (3Ejc), 5; = 7(5, + (3Ejc). (33) 

This mixing of the E and B fields gives the required transformation of the polarization 
directions. In this manner we find R2 in the frame of Fig. [TJ 



rE2 



^2\}CK )W - W2 



rB2 



e2(lCK jW — W2 



IJ,2{icK')w + W2' e2{icK')w + W2 

W2 = a/u'^ + f ^ + €2{iCK') fJ.2i^CK') k'"^ , 
Pe2 = nE2 ® T^E2, Pb2 = ^32 ® 1152, 

/ 

1 

nE2 



nB2 



K^u^ + f ^ — 2ij3Ku — + f^) 



K^v? + f 2 — 2i(3K,u — I3'^{k,'^ + w^) 



f3vw 
-v{k + \l3u) 
y KM — i/3(fi;^ + v'^) j 

I liu^ + f 2) + ^KU \ 

w{u — iISk) 
vw 



(34) 
(35) 

(36) 
(37) 

(38) 
(39) 



Note that the reflection coefficients (35 1 are evaluated using quantities in the frame co-moving 
with plate 2. 



We have verified the solution ([23])-([39]) by comparing it with a brute-force numerical 
solution of the boundary-value problem for the Green tensor. A novel feature of our approach 
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is that the exact analytical solution is obtained purely by physical reasoning, without directly 
solving the differential equation (|7]). 

4. The stress tensor 



We can now compute the stress tensor. It is first necessary to drop the bare Green tensor 
( pO) ) in ( [23] ), as this gives the diverging zero-point stress in the absence of the plates 0. The 
remaining part of the symmetrized Green tensor (16) gives the expectation values ([T7|)- (18), 
which determine the stress tensor ([!]). Some labour is required, as the inverse matrices in (23 ) 



must be evaluated and simplified; the details are given in Appendix A The stress tensor is 



diagonal, so there is no lateral force on the plates. The component cr^x of the stress is constant 
between the plates and is equal to the perpendicular force per unit area F on the plates. A 
positive sign for F means the plates are attracted to each other, whereas a negative sign means 
they repell each other. We introduce the quantities 



-1 -1 2au, 
El' E2^ 

-l»,-l„2at« 
El' B2 



A-EE 
A-EB 

and write the force per unit area 

F = — 

47]-3 



Abb 
Abe 



Bl' B2^ 
Bl' E2 



1, 



X 



du i dvw 

) J — oo 

{Aee + Abb){u'^ + 



iKuPf - {Aeb + Abe)w'^v'^P'^ 



AeeAbb{u^ + v^- l^^u^3Y - AebAbew^v^(3^ 



(40) 



Note that the numerator in (40) becomes the denominator if the sums of pairs of A's are 
replaced by products. Despite the presence of imaginary terms in the integrand, the expression 
( [40] ) is a real number. This is easily seen by considering a Taylor expansion in u of the 
integrand. Since the permittivities and permeabilities are real on the positive imaginary 
frequency axis [l2T]|. this expansion consists of real terms even in u and imaginary terms odd 
in u\ the latter vanish after the integration with respect to u. When /9 = 0, Eqn. (40) reduces 
to Lifshitz's result ill 121 

F{f3 = 0) = —j^ d^j duj dvw{A-^'E + ■ (41) 

The perpendicular force clearly cannot depend on the sign of (3, so the lowest order 
correction to (41) when /3 7^ is proportional to (a lateral force could of course be 



proportional to {3 to lowest order). The expression for this 0{(3 ) correction to (41 ) depends 



on the dispersion in Plate 2, as £2 and /i2 in the reflection coefficients (35 ) are evaluated at 



the co-moving (imaginary) frequency given in (31 ); one must therefore also expand rE2 and 
rB2 to 0(/9^) to obtain the required correction. If we assume there is no dispersion then the 



perpendicular force (40) expanded to 0(/5^) is 

he 

47^3 



/•oo /»oo /»oo 

/ d/t / du I dvw [A^^ + A^^ 

Jo J — 00 J — 00 



+0- 



^^2^2^2aw 



jrEi - rBi){rE2 - ^^2) (e^""' - r eiT E2r biT B2) 
(y. ^2)2 {rEirE2rBirB2AEEABBY 



(42) 
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The Casimir force ( [4T| ) between non-moving plates is attractive {F > 0) for dielectrics 
but can be repulsive {F < 0) if there is a magnetic response. If each plate has either an infinite 



permittivity or an infinity permeability, (40) and (41 1 give the same result so the motion has 
no effect on the force. For general dispersive dielectrics we find that the motion serves to 
add an attractive component to the force compared to the non-moving case, whereas with a 
magnetic response this component can be attractive or repulsive. 

Presentations of the static problem [2| {(3 = 0) do not usually address the issue of the 
electromagnetic stress tensor inside the plates. It turns out that the relevant component of the 
stress tensor (the xx-component) is zero in the plates in the static case so the correct Casimir 
force can be found by restricting the analysis to the stress between the plates. (The yy- and 
2;2;-components of the stress in the material are not zero, however, giving a non-vanishing 
Casimir energy in the plates.) It is important to check whether this situation is maintained in 
the moving case /3 7^ 0. To find the stress tensor in Plate 1 we require the Green tensor for the 



wave equation in the material. The bare Green tensor is a simple generalization of ( 20 ) and the 
physical reasoning used above to write down the Green tensor between the plates can readily 
be applied to find the Green tensor in the material by moving r and r' in Figure 1 into Plate 1 . 



The details are presented in Appendix B and the results show that the electromagnetic stress 



tensor in Plate 1 is diagonal with zero xx-component, as in the static case; hence the complete 



Casimir force on the plate is indeed the perpendicular force (40). By Lorentz invariance the 
same applies to Plate 2. 

As mentioned in the Introduction the previous attempts fil2il at an exact solution of this 
problem failed to calculate the Green tensor correctly (moreover, the important issue of the 
stress tensor in the plates was not addressed). The expressions for the Green tensor and any 
non- vanishing vacuum forces cannot be written as a sum of contributions from two orthogonal 



polarizations as in the static case (41 ). In p2| it was mistakenly concluded that this feature 
of the static result would be preserved in the moving case (this error was also made in ifTOll ). 
As we have demonstrated, the relevant polarization decomposition for a reflecting wave is 



different at each plate, even to first order in /?, due to the Lorentz transformation (32)-(33 1; 
one can only use the same polarization decomposition at each plate if one ignores the motion. 
There are therefore four, rather than two, polarization directions that have to be considered 
for each plane wave. The lack of simple orthogonality relations between these four directions 



is the reason why the correct solution (40) for /? 7^ looks nothing like (41 ). Nevertheless, 



given the complexity of the Green tensor from which it is calculated (see Appendix A), the 



expression (40) has a certain elegance of its own. 



5. Conclusion 

Shear motion at a constant speed of one infinite plate parallel to another modifies the Casimir 
force between them compared to the non-moving case, but it does not induce a lateral force 
on the plates. It has been shown [fT6| that the Casimir-Polder force between a polarizable 
particle and an infinite plate, at rest relative to each other, can be obtained from the Lifshitz 



formula (41 ) for the force between two infinite plates by a limiting procedure wherein one of 
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the plates is allowed to become rarified. Alternatively [|32l . the Green tensor outside a single 
plate can be obtained as a limit of the two-plate Green tensor and used to calculate the force on 
a polarizable particle. It is clear from both these methods that our result implies the absence 
of "quantum friction" on a particle moving a constant speed parallel to an infinite plate. 

Very different is the case of two finite bodies moving past each other with constant 
velocities: this is a complicated dynamical problem, even if the effect of the vaccum forces on 
the velocities of the bodies is neglected. But it is questionable if some feature of the interaction 
of the bodies could be usefully singled out as "friction", especially if the separation of the 
bodies is not small compared to their sizes. In particular, depending on the state of motion, 
a net attractive force between dielectric bodies could act to speed them up rather than slow 
them down. 
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Appendix A. Green tensor between the plates: expanded form 



The Green tensor given by ([23])-([39]) contains the inverse matrices 



(1 - e-^-'^RiRa) ' and (l - e-^-^RaRi) 



-1 



(A.l) 



which must be expanded before computing the stress tensor. These matrices clearly have a 
geometric-series expansion and this, together with the form of the reflection operators Ri and 



R2 (see (24) and (34)) shows that they can be written in a basis constructed from the vectors 



riEi, n^i, n£;2, and 1152 • A little thought leads one to conclude that 

(1 -e 

M2 = C2EE T»-E2 

(1 -e 

CiEE riEi ® rVEl 



-2-"^RiR2) ' = 1 + M2, 



^E2 + C2BB ^B2 ® T^B2 + C2EB T^E2 



nB2 + C2BE nB2 riE2, 



Ml 



'R2R1) 

ClBB n^i 



= 1 + Mi, 

n^l + CiEB ^El 



nsi + CiBE nsi (X) n£;i, 



(A.2) 
(A.3) 
(A.4) 
(A.5) 



where the various c's are unknown coefficients. At this stage (A.2)-(A.5) represent a 
conjecture; it will be verified at the end of the procedure that they are true. The unknown 



c-coefficients in (A.3 ) and (A.5 ) can be found by taking successive dot products of ( A. 1 ) with 
polarization vectors; to give one example, C2eb is given by 



C2EB 



(1 - e-^^-^RiRs) 



-1 



1 



T^B2)j[nE2) 



(A.6) 



We used Mathematica to calculate all of the required dot products of the form (A.6) and 
thereby found all of the unknown c-coefficients. Defining a, A and v by 

e^"'"(rsiri52A^ + rBirB2X^ + rE2rBiJy^ + rEirB2Jy^)] ,(A.7) 



a 



e^""' + rEirE2rBirB2 



A — n^;! ■ nE2 



V 



(A.8) 
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we write the results for the c-coefficients: 

C2EE = OirE2 [-TeiT 31^32 + e^'""(ri5iA^ + tbiv"^)] 



2^ 



C233 = arB2 [-rEirE2r3i + e^""'(rBiA^ + teiI 

CiEE = oirEi [-rE2r3ir32 + Q^""" {r e2>^^ + rB2i 

Ci33 = arsi [-rEirE2r32 + e^''"'(rB2A^ + rE2^'^ 

C23E = ae'^"''"rE2{rEi - ^bi)Az/, C2E3 = ae^'''" r B2{r ei - rBi)Ai/, 

Ci3E = oie^"''"rEi{rE2 - ^52) Ai^, Cieb = ae^'''" r Bi{r e2 - ?^B2)Ai^. 

The expanded form (A. 2 )-( |A.5 ) of the inverse matrices ( |A.l I is now determined and the final 
step in the procedure was to verify using Mathematica that (A.2)-(A.5 1 are indeed correct. 
From (1231) we see that the Green tensor also contains 
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(A.9) 
(A.10) 
(A.ll) 
(A.12) 
(A.13) 
(A.14) 



(1 _ e-2"'"R^R2) and (l - e"2"'"R2Ri) ^Rs- 

Using the above results for ( A.ip , these can be written in the form 

(1 - e-2'""RiR2) Ri = RxNi, 
Ni = diEE n^i ® riEi + diBB n^i ® n^i + diEB ^ei ® n^i + diBE n^i ® riEi, 
(1 - e-2-'^R2Ri)"' R2 = R^Na, 

N2 = d2EE TtlE2 ® Tt\E2 + ^25^ ^^32 ® ^32 + (^2^5 ^32 ® 11^2 + d23E ^32 ® ^E2, 



where R^^ is defined by (25 1 and the (i-coefficients are 

diEE = ae^'^^TEi [e^'^" - rBi{r32>? + rE2i^^ 
diB3 = ae^'^^TBi [e^''"' - rEi{rE2\^ + rB2i^^ 
di3E = diE3 = <ye^"''"rEir3i{rE2 - ^52) Az/, 
d2EE = ae''''"rE2 [e'""" - rB2{r3iX' + VEiiy' 
d233 = ae2"-rB2 [e^'^"' - rE2{rEiX' + tbi^' 
d23E = d2E3 = 0ie^"''"rE2r32{rEi - ^Bi)Az/. 
Finally, we can write the Green tensor ( [23] ) in expanded form 



(A.15) 

(A.16) 
(A.17) 
(A.18) 
(A.19) 

(A.20) 
(A.21) 
(A.22) 
(A.23) 
(A.24) 
(A.25) 



G = Gb + 



+ 



(A.26) 



We also verified directly using Mathematica that ( |A.26p is equal to (23 1. In the form (A.26) 
the Green tensor was used to calculate the electromagnetic stress tensor as described in the 
main text. 



Appendix B. The electromagnetic stress tensor in Plate 1 



In this Appendix we discuss the Green tensor for the vector potential inside Plate 1 and the 
resulting electromagnetic stress tensor. In Plate 1 the Green tensor satisfies (we use a subscript 



No quantum friction between uniformly moving plates 



14 















r 






r 










i 



x=0 x=a 
Figure Bl. Calculation of the Green tensor in Plate 1 . 



m for the Green tensor in the material) 



r,r ,u] 



/ dtG„(r,t,r',0)e'-*, 

^0 



V X 



V X 



UJ 



lS{r 



(B.l) 
(B.2) 



The bare Green tensor in the plate is that for a material filling all space, with no material 
boundaries. In terms of the Fourier transform 



the bare Green tensor is [fT6l 

Grab{x,x',U,V,i^) 



dy 



dz G„(r,r',iOe^ 



-iu{y-y')-iv{z-z') 



q-wi{x-x')q 



X > X 



-jUii (x—x') 



— 5 X X 



1 



i± = 

c 















±Wi 










iu 








iu 








v 


iv 


J 




v 


iv 


J 





Wi 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



The two possibilities in (B.5 ) are plane waves propagating to the right (first line) or to the left 
(second line), with wave vectors ki = (±iwi,u, v). 

The Green tensor Gm(i', ^, r', represents the vector potential at a point r in Plate 1 
resulting from a radiating dipole at another point r' in the plate (see Figure Bl I. Because of 
the Fourier decomposition (B.3 1 we need only consider plane waves. As in Section|3]we can 
simply write down the solution for the Green tensor. Let be the reflection operator for a 
right-moving plane wave in the plate at the a; = boundary and let be the transmission 
operator for that plane wave at the x = boundary. Let Ti be the transmission operator 
for a left-moving plane wave in the vacuum between the plates at the x = boundary. The 
bare Green tensor ( B.5| ) represents direct propagation from x' to x. A plane wave can also 
propagate from x' to x by reflecting off the boundary at x = 0, or by being transmitted through 
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the a; = boundary, reflecting once off Plate 2, then reflecting any even number (including 
zero) of times off both plates before finally being transmitted at the a; = boundary. The 
Green tensor is therefore 

= Gmb + e"'i(^'+^)R„g^+ (B.7) 



mV m-\- • 



(B.8) 

The exponential factors in ( |B.7| ) are deduced from the propagation distances of the waves, as 
described in Section [3]. 

It remains to find the operators R^, T^, and Ti. This is done as for the operators Ri and 
R2 in Section [3]. We require the unit vectors in the directions of the E- and B-polarizations 
in the plate; for the E-polarization the direction is 

/ \ 



1 



U 



y U J 



whereas for the B -polarization we have 

1 



n 



± 

Bm. 



\ 



where n^^ is for a right-moving (left-moving) plane wave. We can now write the 

operators R^, T^, and Ti in terms of scalar reflection and transmission coefficients: 



RjT! 

T 

Ti 



- Em^Em ^ ^E 

- ^Em^Em ® ^Em ~ ^Bm (RajUfil ) 

tElT^El ® riEi + tBin~B^ ® riBl- 



n 



Bmi 



n 



Bmi 



(B.9) 
(B.IO) 
(B.ll) 



The reflection and transmission coefficients in (B.9 )-( B.ll) for the two polarizations are a 
standard result [|27l : 

Wi — HiW Wi — SiW 



rEr, 



tEr. 



^El 



2Wi 



-- 1 - tei, 
1 + Tei, 



— ^Eli r Bm 

Wi + EiW 

2y/e^wi [el 



-rBu 



tBr. 



^Bl 



/ilWl + EifiiW 
2y/€ifliW _ [JI^ 




l-TBl] 



1 + rm), 



Wi + ^iW Wi + EiW 

The Green tensor ( B.7| ) is now completely specified. 

The quantity (1 — e~^"'"R2Ri)~^R2 in (B.7) has been expanded in Appendix A 



(B.12) 
(B.13) 
(B.14) 

so we 



can immediately write the expanded form of the Green tensor: 

Gm = Gmb + e"'^("'+") [R„^ + e-2-"TiR,N2T„] (B.15) 
This is a suitable form for evaluation of the electromagnetic stress tensor. In macroscopic 
electromagnetism the stress tensor is[|] 

^ (B.16) 



£o(D ® E) + i(H ® B) - -t{eo{D ■ E) + /Iq ^(H ■ B)). 



I It is shown in |fT6ll that the Casimir stress in a material is given by (B.I61; this result was challenged in 
but see 1.341 . 
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This is computed from 



(D(r)®E(r)) = -^ / dee(iOe'G^(r, r, i^), (B.17) 



vr 

(H(r)®B(r)) = ^ lim T V x G^(r,r',iOx V', (B.18) 

which follow from ( fTTj ) and ( [TS] ) and the standard definition of the D and H fields. The 
regularization procedure is again to subtract the diverging stress in the absence of the material 
boundary, which in this case means dropping the bare part Gmh of the material Green 
tensor [fT6l . The stress tensor ( |B.16[ ) is found to be diagonal and has zero xx-component, 
as discussed in Section IH 
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